Universal scaling in a strongly interacting Rydberg gas 
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We study a gas of ultracold atoms resonantly driven into a strongly interacting Rydberg state. 
The long distance behavior of the spatially frozen effective pseudospin system is determined by a 
set of dimensionless parameters, and we find that the experimental data exhibits algebraic scaling 
laws for the excitation dynamics and the saturation of Rydberg excitation. Mean field calculations 
as well as numerical simulations provide an excellent agreement with the experimental finding, and 
are evidence for universality in a strongly interacting frozen Rydberg gas. 



I. INTRODUCTION 



The concept of universality appears in many different 
fields of physics [l|, biology Q, economics [|[ and vari- 
ous other systems. It allows to describe the behavior of 
a system without actually knowing all the microscopic 
details of its state. A particular class of universal scaling 
behavior can be found close to second order phase transi- 
tions. The characterization of the corresponding critical 
points in terms of universality classes Q has become cru- 
cial for the understanding of classical as well as quantum 
phase transitions. Quantum degenerate gases can serve 
as a well controlled model system for the exploration of 
universal scaling behavior and quantum phase transitions 
Q in strongly interacting cold atomic systems. Here, we 
show that the experimental data supports the appearance 
of universal scaling in ultracold Rydberg gases, which is 
in agreement with the recently predicted existence of a 
quantum critical point Q . 

The key ingredients of the described experiments are the 
combination of a Rydberg gas in the 'frozen' regime |7[ 
with strong interactions among the Rydberg atoms [8J, 
and the ability to coherently drive the system H as a 
pseudo spin. There exists a variety of interaction mech- 
anisms among Rydberg atoms giving rise to blockade 
phenomena, which are intensively studied [1, [l(| E| ex ~ 
pcrimcntally. Recently several groups also focused on 
the coherent properties of frozen Rydberg gases in the 
regime of weak fl2l. [l3l. [Tij , as well as strong interac- 
tions 0, [TH, Uti Il7l. This unique combination of 
strong interactions with long coherence times led to var- 
ious proposals for quantum information processing using 
Rydberg atoms [Ht M, HI, H3] . 

In this letter, we apply the theoretical framework of 
a quantum critical behavior in strongly interacting Ryd- 
berg gases @ to experimental data (see Fig. QJi), which 
has been previously analyzed with respect to coherent 
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FIG. 1: Universal scaling of the Rydberg fraction in the sat- 
urated excitation regime (a) The saturated Rydberg excita- 
tion is obtained with a laser coupling strength of f2=2-7rxl54 
kHz in dense ultracold atomic clouds with densities n = 
[3.2xl0 19 (< ), 6.6xl0 18 (v), 2.8xl0 18 (n)] m" 3 @j. (b) 
Scanned parameter space for the individual excitation curves 
depicted in the n - Q. plane, (c) Saturated Rydberg fraction 
fn as a function of the dimensionless parameter a — hQ/Cen 2 
for a three dimensional configuration (■) and numerical sim- 
ulations (•). The experimental and numerical data are fit- 
ted (solid lines) to power laws of the form fa ~ a 1//<5 from 
which the critical exponents 1/8 = 0.45 ± 0.01 (exp.) and 
1/5 — 0.404 (num.) are extracted. 



and collective excitation of Rydberg atoms in the strong 
blockade regime @. The relevant parameters of the 
experiment are the density of particles n, the coupling 
strength of the driving laser field f2 (see Fig. Q})), its de- 
tuning from resonance 5l, and the interaction strength 
among the Rydberg states determined in our case by the 
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FIG. 2: Equation of state of a ferromagnet in comparison to 
a strongly interacting Rydberg gas a) Magnetization M of a 
classical ferromagnetic Ising model: The system exhibits a 
critical point at t = 1 — T/T c = and H — 0. The mag- 
netization M is given by a power laws (a) M ~ iP and (B) 
M ~ H 1 / 6 . The exponents have been calculated with the 
Ising model to be /3 = 0.31 (measured: 0.32-0.39) and 5 = 5 
(measured: 4-5) b) Excited state fraction Jr. of a strongly 
interacting Rydberg gas. The universal scaling laws are given 
by power laws of the rescaled detuning A of the driving field 
(?) /i?(a = 0) ~ A 13 and the rescaled coupling strength a as 
depicted by (3) f R (A = 0) ~ a 17 *. Within mean-field the- 
ory the plotted parameter fn is given for a three dimensional 
cloud with van der Waals interaction by a = — A/ 
with 5 = 5/2 and /3 = 1/2. 

van der Waals constant C§. On resonance (Sl = 0), 
these parameters can be merged into a single dimension- 
less parameter a = Ml /Can 2 . We find, that all exper- 
imental data taken from [8| collapse to a simple power 
law as a function of this parameter a, see Fig. [TJ), which 
is in agreement with the predicted universal scaling be- 
havior. This implies that the interpretation in terms of 
a continuous quantum phase transition Q is supported 
by experimental data, and thus provides a firmer theo- 
retical foundation compared to previous scaling models 
Q. While the scaling models essentially depend on two 
parameters (Rabi frequency and ground state density), 
there is only one single parameter in the universal scal- 
ing function, which is the length scale diverging at the 
critical point. We also show that these scaling models 
can be rigorously derived from the mean-field approach 
given in [fj|. 



II. UNIVERSAL SCALING THEORY 

Textbooks on statistical mechanics [l[ often introduce 
universal scaling as a critical phenomenon, which can be 
found near the critical point of a second order phase 
transition. When approaching the critical point the 
system becomes scale invariant. This means that mi- 
croscopic details of the system become irrelevant, and 
the macroscopic behavior is dominated by its long-range 
physics, associated with a diverging length scale £. Math- 
ematically, a function f(s) is called scale- invariant, if 
/(As) = cf(s). Using a series expansion one can see 
that solutions to this equation are given by power laws 
of the form f(s) oc s u . Hence, near the critical point all 
observables can be described by power laws of the diverg- 



ing scale £. Since the critical properties are dominated 
by long-range physics many different systems show the 
same critical behavior. This leads to the classification of 
critical exponents in terms of universality classes, which 
are determined by the spatial dimension, the symmetries 
of the Hamiltonian and the long-range behavior of the 
interactions. 

Phase transitions occur when the free energy of a sys- 
tem shows nonanalytic behavior. In classical systems this 
is always related to a change in temperature. However, 
in quantum systems at T = there is a another possibil- 
ity: the ground state energy can become nonanalytic in 
the case of an avoided crossing with vanishing gap or an 
actual level crossing (23|. Analogously to their classical 
counterparts one can classify quantum phase transitions 
into first, second, and inifinite order transitions. 

A. Universal scaling in a Ferromagnet 

The most prominent example of a second order phase 
transition is a ferromagnet close to the Curie tempera- 
ture Tc- For T > Tc the system is completely demag- 
netized and rotationally invariant. Lowering the tem- 
perature below Tc, the system enters the ferromagentic 
phase characterized by a finite magnetization M , which 
breaks the rotational symmetry. The appearance of the 
magnetization can be described in terms of an order pa- 
rameter. Above Tc the magnetization M is simply zero, 
and exhibits a linear magnetization in the presence of an 
external magnetic field H . The equation of state close to 
the critical point at T c is shown in Fig. [2^) in terms of 
the reduced temperature t = 1 — T /Tc and the external 
magnetic field H. Then the system is described in terms 
of universal scaling laws of the form M ~ t 13 for H = 
and M ~ H 1 / 5 for t = 0, indicated by the curves (a) and 
(fJ), while the diverging length satisfies £ ~ \/t v . 

B. Universal scaling in a strongly interacting 
Rydberg gas 

Now we want to draw the analogy of the magnet close 
to the Curie point to a strongly interacting Rydberg gas 
driven by a laser field. The frozen atomic gas is assumed 
to consist of an ensemble of spatially fixed pseudospins 
with two electronic states, one being the ground state and 
the other a highly excited Rydberg state. The coupling 
CI between the two states is achieved by a monochromatic 
light field with a detuning Sl with respect to the energy 
splitting of the two states. The Rydberg states inter- 
act strongly with a general interaction potential C p /r p , 
which in the present experimental situation is dominated 
by the van der Waals interaction with p = 6. The corre- 
sponding iV-particle Hamiltonian then reads 

i i j<i 
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where the <T x , y ,z are Pauli matrices, the position of 
the atom i and Pee = (1 + Oz 1 ' ) /2 is the projector on the 
excited Rydberg state. It has been recently shown 
that this Hamiltonian features a quantum critical point 
at £1 = 5 = for p > d, where d is the dimensionality of 
the system. In contrast to the above classical example, 
the critical point appears even at zero temperature by 
varying the detuning 5l- Consequently, the role of the 
reduced temperature t is now taken by the dimensionless 
detuning A = fo5i,/E c with the characteristic energy 
E c = C p n p ' d . For A < 0, all atoms remain in the ground 
state for 17 — ► 0, while for A > a finite number Nr of 
atoms are excited into the Rydberg state. The fraction 
of excited Rydberg atoms Jr = Nr/N then plays the 
analog role of the magnetization M in the example of 
a ferromagnetic phase transition discussed above. In 
analogy, we can draw the fraction of excited Rydberg 
atoms fn close to the critical point as a function of the 
detuning A and the parameter a = Ml/E c , see Fig. [5b). 
Note, that the continuous behavior of /r at the critical 
point is a special property of the long-distance behavior 
of the interaction potential. The coherence length £ 
is determined by the characteristic length scale of the 

correlation function ((Pee — fit) (Pee — /k)) 5 which 
gives rise to a diverging length scale close to the critical 
point. This quantity corresponds to the blockade radius 
in the system, i.e., the radius around a Rydberg atom, 
up to which an additional Rydberg excitation is strongly 
suppressed. 



of the full behavior of fa in terms of a general scaling 
function 

fn = a 2d ^ 2p+d ^ y f - ^ (2) 

This is the inverse of the exact mean field result 
a = /il 1 - A /fR P \ shown in Fig. Eb). As the system 
turns classical in the limit a — > 0, the behavior of the 
unknown function x(y) m the limit y — > ±oo can exactly 
be determined to be x(v) ~ U d ^ p and x(y) = 1 / J/ 2 , 
respectively. It is important to note that the mean-field 
result for 5 is the only consistent scaling exponent with 
these limits. This indicates that the classical behavior 
of the system at a — fixes the exponent 5 to its mean 
field value even for low dimensions. In contrast to a 
classical critical point, for a quantum critical point the 
dynamical behavior is coupled to the static properties 
via the dynamical critical exponent z, i.e. r ~ £ z where 
t describes the characteristic time scale close to the 
critical point. Here, we find the dynamical critical 
exponent to be z = p, which implies that the relaxation 
is domitated by the frequency \A/V&Q in agreement with 
the above superatom picture. 



III. EXPERIMENT 
A. Universal scaling in an inhomogeneous sample 



It remains to identify the critical exponents f3 and 
1/5 for the universal scaling of the observable parame- 
ter f R ~ and f R ~ a 1/s . For a = the Hamilto- 
nian |T]) is classical and by minimizing its energy one 
obtains (3 = d/p. The correlation length £ is deter- 
mined by the averaged spacing between the Rydberg 
atoms via £ ~ a/f^ d (here, a^ 1 = tfn denotes the 
averaged intcrparticlc distance of a c?-dimcnsional sys- 
tem). To obtain a value for 1/6, it is useful to take a 
closer look at the excitation dynamics of a strongly in- 
teracting Rydberg system (a -C 1) driven at resonance 
(A = 0), which also complies with the experimental sit- 
uation When a approaches zero, the system enters 
the strongly blocked regime with Jr <C 1. Within the 
blockade radius £ only one excitation is shared over a 
large number of atoms N), ~ n£ d resulting in a collec- 
tive state \ip e ) = ^Eil b i 151,52, 9z, ...,e u ...g Nb ), which 
shows an accelerated temporal evolution with a collec- 
tive Rabi frequency y/N^n @, [H, H3- In this so called 
superatom model the blockade radius can be evaluated 
by equating the interaction energy with the collective 
coupling strength as C p /£_ p = \ZT7bQ. Note, that £ di- 
verges for ct — ► 0. This procedure determines 1/5 to be 
1/6 = 2d / (2p+ d) and the case for d = 3 and p = 6 is car- 
ried out in some more detail in [2f| . The same result for 
1/5 is obtained using standard mean-field theory 0]. In 
addition, the mean-field solution leads to a description 



In the actual experiment, the atoms are well de- 
scribed as a thermal gas trapped by a harmonic potential. 
Then, the three dimensional density distribution of the N 
ground state atoms has a Gaussian shape with radii given 
by the standard deviations (J x ,y,z = y/ksT '/2mu Ii!( . 2 , 
which are determind by the trapping frequencies uJ x ,y,z, 
the mass m, the temperature T of the cloud and Boltz- 
manns constant kg 

(2TxYl 1 a x a y a z \ 2o% 2af j 2af ) 

Note, that the temperature T associated with the ki- 
netic energy of the atoms is decoupled from the dynam- 
ics of the Rydberg excitations in the frozen Rydberg 
gas. Within the local density approximation, we can de- 
scribe the properties of the system by a local parameter 
a(r) = Ml/ Cqti(y) 2 and the total Rydberg fraction Jr, is 
given by (d = 3, p = 6) 

f R = — [ dr 3 fR(r)n(r) - -i — -^a 1 ' 8 . (4) 

JR N J (-2/ 5+ If' 2 V ; 

Here, a is the peak value in the trap center a = 
M7 / Cqii(0) 2 . Consequently, we find that the critical ex- 
ponent 5 is not modified by the harmonic trapping poten- 
tial within the local density approximation, and reduces 
to the value given in the thermodynamic limit. 
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FIG. 3: Universal scaling behavior of the excitation rate. The 
rescaled excitation rate ga for a three dimensional density 
distribution is shown for experimental data (■) and the cor- 
responding numerical simulation (•)). A linear fit to a power 
law gn ~ a 1 results in a critical exponent of 7 = 1.25 ± 0.03 
(exp) and 1.15 (mini). 



TABLE I: Comparison of the different results for the critical 
exponents 7 and 1/5 for a three dimensional as well a one 
dimensional density distribution. The theoretical values are 
given in Eq. (fSJ) and Eq. The results of the numerical 

simulation have been achieved by integrating the Hamiltonian 
(UJ for up to 100 particles [(|. The experimental results are 
obtained by fitting power laws to the data shown in Fig. [T] 
and Fig. 





7 (g R ~ a 1 ) 


1/6 (fR-a 1 ' 6 ) 


experiment [Id] 


1.08 ±0.01 


0.16 ±0.01 


theory 


14/13 » 1.08 


2/13 » 0.15 


numerical simulation 


1.06 


0.150 [6] 


experiment [3d] 


1.25 ±0.03 


0.45 ±0.01 


theory 


6/5 = 1.2 


2/5 = 0.4 


numerical simulation 


1.15 


0.404 [6] 



C. Experimental results 



B. Experimental setup and procedure 



A detailed des crip tion of the experimental setup can 
be found in [H, l26l | and here only a rough outline of 
the experimental procedure is given. First we prepare 
a magnetically trapped cloud of Rubidium atoms spin 
polarized in the 55*1/2 state. The atomic cloud has a 
temperature of 3.4 /xK and peak densities no close to 
10 20 m~ 3 . This corresponds to phase-space densities 
below quantum degeneracy to avoid a bimodal density 
distribution [25j of a Bosc-Einstcin condensate. To alter 
the atomic density without changing the size of the 
cloud we use a Landau-Zener sweep technique Q to 
alter the number of trapped atoms from N = 1.5 x 10 7 
to N — 5 x 10 5 . In this process the temperature remains 
unchanged and consequently the physical dimensions of 
the cloud. With the known harmonic oscillator potential 
of the Ioffe-Pritchard type trap, all parameters of the 
atomic clouds are known. The excitation to the 43S , 1 / 2 
Rydberg state is done with a resonant two-photon 
transition via the 5P3/2 state. To avoid population of 



the intermediate 5P 3 / 2 state the light is detuned with 
respect to this state by 6l=2ttx478 MHz to the blue. 
The coupling strength £1 of this effective two-level system 
is altered from 27rx31 kHz to 27rxl54 kHz by adjusting 
the laser intensity. The excitation dynamics are investi- 
gated by a variation of the excitation time from 100 ns 
to 20 /is, which is short compared to the excited state 
lifetime of 100 /is. After excitation the Rydberg atoms 
are field ionized and the emerging ions are detected 
with the help of a multi-channel plate. The resulting 
strongly blockaded excitation dynamics is determined by 
the strong repulsive interaction among Rydberg states, 
which is in our case given by the isotropic van der Waals 



interaction with Ce, 



1.7 x 10 19 a.u. for the 43S state. 



In the experiment a is changed non-adiabatically by 
switching on abruptly the coupling laser field f2. Then, 
the number of excited Rydberg atoms Nft(t) undergoes 
a dynamical evolution, which saturates in Nr as shown 
in Fig. []Ji). The inital increase in the number of Ryd- 
berg atoms is well described by a rate R, and this relax- 
ation time is experimentally deduced by a fit of the time- 
evolution of the Rydberg excitation Nn(i) by an expo- 
nential saturation function Nr(£) = Nr(1 — e~ Rt ' NR ), 
which allows to extract both the initial excitation rate 
R and the saturation level Nr. In the strongly blocked 
regime (a <C 1) the rate R is determined by the collective 
Rabi frequency and the saturation level Nr, which 

is close to that of the ground state of the system. Previ- 
ously, we have examined the data and its dependence on 
the ground state density n and the Rabi frequency O us- 



ing a general expression R ~ n KR Q R and Nr 



[8|. The scaling behavior for a variation of the ground 
state density n and the coupling strength CI gave a strong 
evidence for a coherent collective excitation dynamics in 
the strong blockade regime. 

In the following, we analyze these results in terms of 
a universal scaling behavior at resonance (A = 0) by 
rescaling the measured quantities to a dimensionless rate 
gR and and a dimensionless saturation level /r as 



HR 



9r 
Jr 



NC p nP/ d 



7 



Nr 
N 



a 



1/8 



2(p + d) 
2p + d ' 
2p + d 
2d ' 



(5) 
(6) 



First, we would like to point out, that the data collapse 
in to a algebraic relations as shown Fig. QJ) and Fig. 
[3]) is in agreement with the predicted scaling laws. It is 
worth to mention that the numerical simulations based 
on only 10 2 pscudospins scales up to the experimental 
situation with atom numbers of up to 10 7 . For a 
more quantitative analysis of the experimental data, 
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it is important to point out that the radial Gaussian 
radius of the cigar shaped cloud is a x y = 8.6 [im. This 
width is comparable to the blockade radius between two 
Rydberg atoms of roughly 5 /im and places by this the 
geometry in a crossover regime between one dimension 
and three dimensions. Therefore we analyze the data 
wether the experimental setup is better described in 
one dimension with a line density (a = HQ/C^n 6 ) 
or in three dimensions (a = MI/Cqtl 2 ). Fitting the 
observed power laws of the form In ~ 7 In a + c g 
and ln/^j ~ 1/5 In a + Cf we extracted the individual 
exponents. The results are summarized in table |T] and 
compared to the expected theoretical values as well as to 
numerical simulations; the procedure for the numerical 
analysis is described in detail in Q ■ The relatively small 
error bars of the fitted exponents 1/5 show the excellent 
agreement with power laws over a large range in a. 
Nevertheless, the dimensionality can not unambiguously 
assigned, although there is a better agreement for the 
one dimensional case. Here the exponent is dominated 
by the variation in density, which scales to the 6 th 
power. In future experiments it might be possible to 
create better defined dimensionality by adjusting the 
shape of the atomic cloud and/or the strength of the 
interaction by choosing adequate Rydberg states. 



ing theories. This result confirms that the description 
with the effective spin Hamiltonian given in this article 
is correct to a large extent. Therefore the observation of 
a quantum critical point in this system should be within 
experimental reach. This could be done by measuring the 
excited Rydberg fraction when approaching the quantum 
critical point adiabatically starting from a non-critical re- 
gion in the a — A parameter space. Another way would 
be the observation of a crystalline correlation function 
of the excited Rydberg atoms by either a spatial depen- 
dent observation of the Rydberg atoms or in the Fourier 
space by a Laue diffraction experiment with a four wave 
mixing technique (2?| . The measurement of the critical 
exponent is not in complete accordance with a one dimen- 
sional or three dimensional situation, which is most likely 
due to finite size effects. The general form of the scaling 
exponents allows also to apply the simple model to dipo- 
lar systems, which are widely realized in frozen Rydberg 
gases and are also feasible in the context of ultra-cold 
dipolar molecules [28| . 
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